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ABSTRACT 

We use data from the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey to study how 
the spatial variation in the stellar populations of galaxies relate to the formation of galaxies at 1.5 < 
z < 3.5. We use the Internal Color Dispersion (ICD), measured between the rest-frame UV and optical 
bands, which is sensitive to age (and dust attenuation) variations in stellar populations. The ICD 
shows a relation with the stellar masses and morphologies of the galaxies. Galaxies with the largest 
variation in their stellar populations as evidenced by high ICD have disk-dominated morphologies 
(with Sersic indexes < 2) and stellar masses between 10 < Log M/M 0 < 11. There is a marked 
decrease in the ICD as the stellar mass and/or the Sersic index increases. By studying the relations 
between the ICD and other galaxy properties including sizes, total colors, star-formation rate, and 
dust attenuation, we conclude that the largest variations in stellar populations occur in galaxies where 
the light from newly, high star-forming clumps contrasts older stellar disk populations. This phase 
reaches a peak for galaxies only with a specific stellar mass range, 10 < Log M/M 0 <11, and prior to 
the formation of a substantial bulge/spheroid. In contrast, galaxies at higher or lower stellar masses, 
and/or higher Sersic index (n > 2) show reduced ICD values, implying a greater homogeneity of 
their stellar populations. This indicates that if a galaxy is to have both a quiescent bulge along 
with a star forming disk, typical of Hubble Sequence galaxies, this is most common for stellar masses 
10 < Log M/Mq <11 and when the bulge component remains relatively small (n < 2). 

Subject headings: galaxies: evolution — galaxies: stellar content — galaxies structure — galaxies: 
morphology 


1. INTRODUCTION 

In the local universe, galaxies separate into two broad 
classes (Kauffmann et al. 2003; Baldry et al. 2004): (1) 
large disk-dominated galaxies show ongoing star forma¬ 
tion rates (SFRs) comparable to their past averages as 
evidenced by their overall blue colors, and frequent spi¬ 
ral arms, and (2) spheroidal, or bulge-dominated galax¬ 
ies with current SFRs much less than their past aver¬ 
ages. However, especially in the more distant universe, 
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there exists a third class of galaxies which exhibits a large 
amount of stellar population diversity and irregular mor¬ 
phologies (e.g., Giavalisco et al. 1996; Lowenthal et al. 
1997; Dickinson 2000; Papovich et al. 2005). The mech¬ 
anism by which a galaxy transitions from a uniformly 
blue, star forming galaxy (with either disk-dominated or 
irregular morphologies) into a bulge-dominated, quies¬ 
cent galaxy remains an open question. 

Many studies have added to a growing body of evidence 
that the irregular morphologies are due, at least in part, 
to heterogeneous (clumpy) star formation (e.g., Adamo 
et al. 2013; Elmegreen et al. 2005, 2008, 2009a,b; Guo 
et al. 2012; Wuyts et al. 2011a, 2012, 2013; Guo et al. 
2014). Believed to be formed through gravitational in¬ 
stabilities in gas-rich disks (e.g., Keres et al. 2005; Dekel 
& Birnboim 2006; Dekel et al. 2009b; Dekel & Burk- 
ert 2014), these clumps typically contribute 10 — 20% of 
the total galaxy light (Forster Schreiber et al. 2011) and 
most have bluer colors and elevated specific star forma¬ 
tion rates (sSFR) compared to the surrounding regions 
(Guo et al. 2012, 2014). This contrast between a young 
star forming region and the older underlying galaxy pop¬ 
ulation may lead to a large amount of variation in the 
internal colors of a galaxy. 

The ability to investigate the separation (heterogene¬ 
ity) of the different stellar populations (blue star forming 
clumps versus the older, redder stellar disk) becomes in¬ 
creasingly difficult at high redshift, because a galaxy’s 
morphology often depends heavily on the bandpass in 
which the galaxy is observed. For example, observed 
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optical light (rest-frame UV) traces the distribution of 
unobscured star-forming regions (Dickinson 2000) which 
are often non-uniformly distributed. 

Recent observations of clumpy galaxies (e.g., Guo et al. 
2012 ; Wuyts et al. 2012 , 2013; Adamo et al. 2013; Guo 
et al. 2014) find that clumps located toward galaxy cen¬ 
ters are older and redder than clumps located toward 
the outskirts. This suggests that highly star forming 
clumps may migrate toward the center where they could 
condense to form a psuedo-bulge (e.g., Ceverino et al. 
2010). The formation of such a bulge, in turn, could be 
sufficient to stabilize the stellar disk against new star- 
formation (Martig et al. 2009), homogenizing the stellar 
populations and reducing the variations in internal color. 
If this scenario is correct, we would expect to observe a 
relation between the internal color dispersion of a galaxy 
and the formation of a bulge. 

To investigate any potential relation, we require a tool 
to quantify the stellar population diversity of a galaxy. 
Traditional morphological classification methods: visual 
inspection, Sersic profile fitting, concentration and asym¬ 
metry values (Conselice et al. 2004), and Gini coefficients 
(Lotz et al. 2004) are calculated in a single bandpass and, 
as such, do not provide enough insight into the stellar 
population variations. Papovich et al. (2003) developed 
a differential morphological measure, the Internal Color 
Dispersion (ICD), which quantifies the heterogeneity of 
stellar populations (such as with clumps) between band- 
passes; galaxies with non-uniform distributions of young 
stellar populations give rise to a larger ICD than a galaxy 
with a homogeneous or evenly mixed stellar population. 
Similarly, galaxies with variable amounts of extinction 
will also exhibit large ICD values, due to the spatial vari¬ 
ations in the reddening. Combining the ICD with other 
morphological indicators that quantify the distribution 
of light in a galaxy (e.g. Sersic, Gini-M 2 o, concentration 
and asymmetry) provides insight into galaxies physical 
properties such as clumpy star-formation and their evo¬ 
lution in both overall color and structure. 

In this paper, we build upon previous studies (e.g., 
Papovich et al. 2003, 2005; Bond et al. 2011; Law et al. 
2012 ) and we use the ICD as the primary metric in a 
comparative study of the rest-frame UV and optical mor¬ 
phological properties of the largest sample yet (> 800) 
ofl. 5 < 2 :< 3.5 galaxies to date. In addition to broader 
IR coverage, we combine optical and near-IR data from 
HST from the Southern region of The Great Observa¬ 
tories Origins Deep Survey and Cosmic Assembly Near- 
Infrared Deep Extragalactic Survey, which also includes 
the deeper data in the Hubble Ultra Deep Field. We 
investigate the presence of heterogeneous stellar popula¬ 
tions in a sample of 1.5 < z < 3.5 galaxies and inter¬ 
pret our results in the context of galaxy evolution where 
a galaxy transitions from a disk-dominated to a more 
bulge-dominated morphology while experiencing a sup¬ 
pression of the diversity in its stellar populations. 

In Section 2 we discuss the data and sample selection. 
We discuss the ICD measuring technique in detail in Sec¬ 
tion 3. We compare ICD and other galaxy properties in 
Section 4, and discuss possible causes of high ICD in 
Section 5. We discuss the role of star-forming clumps in 
galaxy evolution in Section 6 , and finally summarize our 
results in Section 7. 

Throughout this paper, we use a concordance cosmo¬ 


logical model (Ua = 0.7, = 0.3, and Hq = 70 

km s - 1 Mpc -1 ), assume a Chabrier initial mass func¬ 
tion (IMF; Chabrier 2003), and use AB magnitudes (Oke 
1974). 

2 . DATA 

For this work, we make use of deep imaging data taken 
as part of the Cosmic Assembly Near-Infrared Deep Ex¬ 
tragalactic Survey (CANDELS; Pis Faber and Ferguson; 
Koekemoer et al. 2011; Grogin et al. 2011) of the South¬ 
ern region of The Great Observatories Origins Deep Sur¬ 
vey (GOODS, hereafter referred to as GOODS-S, Gi- 
avalisco et al. 2004) along with extremely deep observa¬ 
tions of the Hubble Ultra Deep Field (HUDF). The large 
variation in depth between the GOODS-S observations 
and the HUDF ensures that we have access to galaxies 
with a large range in stellar mass. 

We take advantage of the third version of the Advanced 
Camera for Surveys (ACS) mosaics in GOODS-S which 
includes images taken with F435W, F606W, F775W, and 
F850LP (hereafter referred to as B 435 , V 606 , ^ 775 , and 
^850 )• la addition, we use imaging obtained from Wide 
Field Camera 3 (WFC3) in F098M, F105W, F125W, and 
F160W (hereafter as Yogs, Y 105 , J125, and Hi 60 ). The 
WFC3 imaging is split into two parts; the deep portion 
covers the central 50% of GOODS-S with the wide por¬ 
tion covering the Southern-most 25% to a depth approxi¬ 
mately 1 magnitude shallower than the deep portion. To 
complete our WFC3 coverage of GOODS-S, we make use 
of the WFC3 Science Organizing Committee’s Early Re¬ 
lease Science program (PI O’Connell; Windhorst et al. 
2011 ) which covers the Northern 25% in Yogs, J125, and 
Hi 60 - 

In addition, we include very deep imaging of the HUDF 
taken with both ACS (Beckwith et al. 2006) and WFC3 
as part of HUDF09 (PI Illingsworth; e.g., Bouwens et al. 
2010; Oesch et al. 2010) and UDF 12 (PI Ellis; Ellis et al. 
2013). Both the ACS and WFC3 datasets contain imag¬ 
ing using the previously mentioned seven bandpasses. 

2 . 1 . Source Detection and Photometric Data 

We use a multi-wavelength photometric catalog of the 
GOODS-S region provided by Guo et al. (2013). In 
this catalog, the sources are detected in a H 160 mosaic 
drizzled to 0.06'ypixel. Along with the source detec¬ 
tions, photometry covering 0.4-8 /am (VIMOS U, B 435 , 
V 606 , 2775 , ^850, Yogs, Y 105 , Jl25, Hl 60 , ISSAC Ks, 
IRAC 3.6/im, 4.5/im, 5.8/im, and 8/im) is provided. The 
photometric data for the HST bands are computed by 
SExtractor (Bertin & Arnouts 1996) on images which 
have had their point spread function (PSF) convolved 
to match the Hieo resolution, while the lower resolution 
ground based and Spitzer /IRAC imaging is processed by 
TFIT (Laidler et al. 2006). See Guo et al. (2013) for full 
details on the construction of the photometric catalog. 

2 . 2 . Photometric Redshifts 

We make use of the photometric redshifts for the galax¬ 
ies in our sample derived by Dahlen et al. (2013). This 
photometric redshift catalog is constructed by combining 
redshifts derived by 11 different investigators using differ¬ 
ent photometric redshift codes e.g., EAZY (Brammer et al. 
2008), HyperZ (Bolzonella et al. 2000), LePhare (Arnouts 
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et al. 1999; Ilbert et al. 2006), etc. These photometric 
redshift codes derive redshifts by fitting different combi¬ 
nation of spectral energy distribution templates and pri¬ 
ors including available spectroscopic redshifts. We make 
use of spectroscopic redshifts where available (see Dahlen 
et al. 2013 for the compilation of spectroscopic datasets). 
The scatter between the photometric redshifts and cor¬ 
responding spectroscopic redshifts, which span ^ ~ 0 — 6, 
is rms[Az/(l + z spec )} = 0.03 where Az = z spec - z phot . 

2.3. Sample Selection 

Papovich et al. (2003) showed that galaxies with het¬ 
erogeneous stellar populations have high ICD between 
their rest-frame UV and optical light. For our galaxy 
sample with 1.5 < z < 3.5, the CANDELS ACS and 
WFC3 bands probe rest-frame UV and optical wave¬ 
lengths respectively. Using the photometric redshifts de¬ 
scribed previously (Section 2.2) we select all galaxies with 
photometric or spectroscopic redshifts in this range. We 
apply a magnitude limit of 25 mag in Hi 60 , which al¬ 
lows for a robust determination of morphological param¬ 
eters and the ICD values (see Section 3.2 for reasoning). 
To avoid stars, we remove objects with stellarity index 
> 0.78; stellarity derived using SExtractor gives a like¬ 
lihood of an object being a star. Applying this set of 
criteria yields a parent sample 3369 galaxies. 

2.4. Stellar Masses 

Because the photometric redshifts are the result of the 
combination of of the probability density functions from 
many routines, Dahlen et al. (2013) show they provide 
more accurate estimates of the redshift. The method 
employed by Dahlen et al. (2013) does not provide ad¬ 
ditional physical properties for the CANDELS galaxies. 
We derive stellar population parameters for the galaxies 
in our sample using FAST (Kriek et al. 2009), a stellar 
population synthesis modeling and fitting code. We fix 
the galaxy’s redshift at the photometric redshift from the 
Dahlen et al. (2013) catalog and fit 14-band galaxy pho¬ 
tometry covering 0.4 — 8 /am with model spectral energy 
distributions to estimate the stellar masses (along with 
dust attenuation) for the galaxies in the sample. We 
used models for a range of stellar population properties 
from the 2003 version of the Bruzual & Chariot models 
(Bruzual & Chariot 2003) stellar population synthesis 
models allowing for a range of attenuation (A v = 0 — 4 
mag) using the Calzetti et al. (2000) extinction law. We 
opt to use Bruzual & Chariot (2003) to facilitate the com¬ 
parison to other studies and because recent work shows 
that the Bruzual & Chariot 2007 models may overesti¬ 
mate the contributions of late type giants to the near-IR 
(e.g., Zibetti et al. 2012). We assumed models with a de¬ 
layed, exponentially declining (Log r = 7 — 11) star for¬ 
mation history, solar metallicity and a Chabrier (2003) 
IMF. 

Using a Salpeter IMF (Salpeter 1955) would to first or¬ 
der increase systematically the stellar masses by ~ 0.25 
dex. Using different assumptions for the stellar popu¬ 
lation metallicities will affect the derived stellar masses 
by ~ 0.2 dex (Papovich et al. 2001). The average un¬ 
certainty in stellar mass for our sample is 0.11 dex, al¬ 
though errors from systematics in the metallicities, star- 
formation histories and other aspects of the spectral en¬ 


ergy distribution fitting likely dominate (e.g., Kriek et al. 
2009; Papovich et al. 2014). 

The errors on the photometric redshifts do not influ¬ 
ence greatly our conclusions. While, photometric error 
uncertainties can affect the derivation of stellar popula¬ 
tion parameters, for small photometric redshift errors, 
like those for our sample, the contribution to the stellar 
mass uncertainty is minimal because small changes in 
redshift translate to small changes in distance (see dis¬ 
cussion in Dickinson et al. 2003). For example, our uncer¬ 
tainty of rms[A z /(1 + z)\ = 0.03 corresponds to a change 
in luminosity distance of ^5% for our sample. Because 
the stellar mass scales quadratically with luminosity dis¬ 
tance, the mass uncertainty from photometric redshift er¬ 
rors is similarly small, ^10%. This error is much smaller 
than systematic uncertainties in stellar mass arising from 
model fitting (see discussion in Papovich et al. (2014)), 
and so we neglect errors from the photometric redshifts 
on the stellar masses. 

Along with stellar masses FAST provides an estimate of 
the total attenuation (Ay) due to dust in a galaxy. We 
use this value in Section 5.5 to understand better how 
the ICD might be effected by reddening due to dust. 

2.5. Morphological Indicators 

We use the structural catalog provided by van der Wei 
et al. (2012) which provides Sersic index, n, (Sersic 1963) 
and the effective radius, R e (measured along the major 
axis of the galaxy), for all sources in the GOODS-S cat¬ 
alog from Guo et al. (2013). These parameters are pro¬ 
duced by GALFIT (Peng et al. 2002) using the best-fitting 
Sersic models for the objects in the Hi6o imaging. 

We also use the non-parametric indicators Gini param¬ 
eter (G) and the second-order moment of the brightest 
20% of the galaxy light, M 20 . These indicators quan¬ 
tify the spatial distribution of light without assuming a 
specific functional form when classifying the morpholo¬ 
gies of disturbed or irregular galaxies. For a complete 
description of this method see Lotz et al. (2004). 

2.6. Color Gradients 

We use color gradient measurements from Hilton et al. 
(in prep), where Hilton et al. compute the color gradi¬ 
ents between the core and outskirts by centering annular 
apertures on each galaxy. They define concentric, ellip¬ 
tical, annular apertures such that the total H 160 light 
of each galaxy is divided equally between three annuli. 
The galaxy core corresponds to the “inner” aperture and 
the galaxy outskirts correspond to the “outer” aperture. 
The color gradient of the galaxy is denoted by 

A(mi - rn 2 ) = (mi - m 2 ) outer - (mi - rn 2 )inner (1) 

where, m\ and m 2 are the respective apparent magni¬ 
tudes. A galaxy where the core has a redder color than 
the outskirts will have a negative color gradient and a 
galaxy with a bluer core than outskirts will have a posi¬ 
tive gradient. 

2.7. Star Formation Rates 

In this work, we measure the total star-formation rate 
(SFR) as the sum of the SFR implied by the observed 
UV luminosity (uncorrected for dust extinction, SFRuv) 
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and by the IR (as estimated from the 24/rm luminos¬ 
ity, SFRir). We convert the observed 24/im luminosity 
into a SFR using the empirical relationship from Wuyts 
et al. (2011b). For galaxies which are not detected in 
24/im, the SFRs are based on an attenuation-corrected 
UV flux using attenuation value estimated by FAST, so 
SFRxotai = SFRuv x Correction. In both instances, 
we follow Kennicutt (1998) to convert the luminosity 
to SFR. The Kennicutt (1998) conversion is Salpeter- 
based, whereas the masses and other derived quanti¬ 
ties are Chabrier-based. To account for the difference 
in IMF, we divide the total SFR by 1.66, SFR Tota i = 
(SFRuv+ SFRir)/1.66 

2.8. Discrete Star Forming Regions 

In Section 5.6 we compare the internal color disper¬ 
sion measurements to the frequency of discrete star form¬ 
ing regions (colloquially “clumps”) in galaxies using the 
clump catalog of Guo et al. (2014). In order to facili¬ 
tate clump finding, Guo et al. (2012) require galaxies to 
be approximately face-on (b/a > 0.5), reasonably ex¬ 
tended (R e > 0.2"), and with a sufficient brightness 
(Hi 60 < 24.5 mag). 

Individual clumps are selected from the UV luminos¬ 
ity contrast between the individual clump and the to¬ 
tal galaxy. From the catalog, we select clumps with 
UV luminosity > 1% of the total galaxy UV light 
(L. dump / L g a i aX y > 0.01) contained within galaxies in our 
sample. While Guo et al. (2014) selects clumps with 
UV luminosity > 8%; we lower this threshold to include 
more UV-faint clumps as they are redder and provide 
a more comprehensive view of the ICD. This provides 
clump measurements for 180 objects, containing between 
zero and eight clumps per galaxy. 

3. THE INTERNAL COLOR DISPERSION 
3.1. Definition 

First developed by Papovich et al. (2003), the Internal 
Color Dispersion (ICD) is a flux-independent statistic, 
defined as the ratio of the squared difference of the image 
flux-intensity values about the mean galaxy color to the 
square of the total image flux, or, 

Npix N P ix 

£ (h,i - a Ii,i -P?~ £ (B 2 ,i - a B m ) 2 

- 

iV pzx iy pzx 

£ {h,i — P) 2 — £ (Bj.i-aBj.i) 2 

i=l i= 1 

x 100%. (2) 

The ICD assesses the morphological difference of a galaxy 
between two passbands. The rest-frame UV-optical ICD 
quantifies the spatial homogeneity (or lack thereof) be¬ 
tween young stellar populations (dominating the rest- 
frame UV light) and older stellar populations, contribut¬ 
ing most of the rest-frame optical light. Galaxies with 
spatially segregated stellar populations with different col¬ 
ors would have large ICD (see Papovich et al. 2003). 
Whereas, a galaxy with a homogeneous stellar popula¬ 
tion throughout would have a low ICD. 

The terms in Equation 2 sum over i = 1 ...N P i X with 
N pix being the number of pixels associated with each 
galaxy, I\ and I 2 are the image flux-intensity values for 


each object obtained in each of two passbands and B\ 
and P >2 are the pixel intensity values from a contiguous, 
similarly sized, region of the two passbands not associ¬ 
ated with any object. 

To determine which pixels belong to a galaxy we de¬ 
fine the aperture size as r p (r] = 0.2), where r p is the Pet¬ 
rosian radius of I 2 as defined for the Sloan Digital Sky 
Survey (Blanton et al. 2001), and rj(r) is defined as the 
ratio of the galaxy surface brightness, /(r), averaged over 
an annulus of radius r, to the mean surface brightness 
within this radius, (/(r)). The Petrosian radius depends 
on the surface-brightness distribution of the galaxy and 
is thus independent of redshift or systematics in the im¬ 
age calibration. The chosen aperture encompasses most 
of the light from the galaxy without adding excessive 
background light. Therefore, the summations are over 
the individual pixels associated with the object or blank 
region of a similar size. 

The scaling factor a is the ratio of total fluxes between 
images I 2 and U, while the linear offset /3 adjusts (if nec¬ 
essary) for differences in the relative background levels 
of the two images. In practice, we calculate these terms 
by minimizing the statistic, y 2 = [{1 2 — a I\ — /?)/cr] 2 , 
adjusting a and f3 as free parameters (where a represents 
the uncertainties on the I 2 flux-intensity values). 

The treatment of the background terms (B 1 and B 2 ) 
are important to the ICD, as the background itself can 
contain an intrinsic amount of internal color dispersion 
based on our definition. Therefore, to eliminate ICD con¬ 
tribution from the background, B 1 and B 2 are selected 
for each individual galaxies to contain the same number 
of pixels (N p i x ) at random locations excluding any preex¬ 
isting objects, and include the same pixels in both image 
1 and 2. In order to minimize the ICD variations due to 
background, we select the median ICD value computed 
from using nine different background regions. However, 
simple fluctuations may not be the only source of ICD. 
Because, the WFC3 imaging PSF has been convolved 
(see Section 2) to match that of the ACS, there exists 
the possibility of correlated noise in the background mea¬ 
surements. This could have an effect if the background 
region is not selected in the same way as the galaxy (e.g., 
choosing individual pixels at random versus a single con¬ 
tiguous region). We test for possible effects by choosing 
individual pixels not associated with an object at random 
from the imaging, choosing a single contiguous square 
region and contiguous circular region. We find no effect 
based on the morphology of chosen background pixels. 

3.2. ICD Dependence on Signal-to-Noise 

Previous studies (Papovich et al. 2003, 2005; Bond 
et al. 2011) noted that the ICD produces unreliable re¬ 
sults for galaxies with low signal-to-noise (S/N). We test 
the behavior of the ICD using simulated galaxies with a 
priori known (fixed) total flux and ICD. Each simulated 
galaxy is placed throughout the HST images at a ran¬ 
dom location and is then checked to insure no overlap 
with an actual object by comparing the location with 
the SExtractor segmentation image. The flux of the 
galaxy is slowly decreased to a level barely above the 
background. We then derive an ICD value for each sim¬ 
ulated object using the definition in Equation 2. 

Figure 1 shows the fractional change (true minus ob¬ 
served divided by true) in the ICD as a function of S/N = 
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Fig. 1.— The fractional change of the Internal Color Dispersion 
as a function of signal-to-noise for two different simulated objects. 
Gray points show individual S/N measurements. Red dashed lines 
show the median value. Green shaded regions shows the area where 
the fractional change is less than 33% from the true value and 
serves as an aid to gauge accuracy. For the low ICD galaxy (top 
panel) we find that galaxies with S/N < 10 in the bluer band the 
median ICD deviates more than 34% away from the true value. 
Therefore, we remove any object from our initial sample that have 
S/N > 10 in the bluer band. The S/N requirement for galaxies 
with intrinsically high ICD could be set to a lower threshold, but 
we choose to constant criterion to keep the selection uniform. 


f / (j2 ’> where / represents the flux-intensity val¬ 

ues, and a 2 is the image variance due to the noise of 
the background. The sum is taken over all pixels iden¬ 
tified as belonging to the object, by computing the Pet¬ 
rosian radius and comparing the enclosed pixels to the 
SExtractor segmentation image. The four simulated ob¬ 
jects illustrate how different levels of ICD ( 2 %, 6 %, 13%, 
and 26%) are effected by the galaxy’s S/N. 

The behavior of the fractional change is dramatically 
different depending on whether a galaxy has intrinsically 
“low” or “high” ICD. For all galaxies with very high S/N 
we recover the expected ICD to within 0.1%. However, 
as the S/N decreases Equation 2 systematically underes¬ 
timates the ICD. This bias is most apparent for galaxies 
with “low” intrinsic ICD. 

Because of the systematic underestimation of the ICD 
it is important to clip the sample at a S/N which still 
provides accurate results. We choose S/N = 10 such that 
the bias is does not impact our results. For example, the 
bias is the fractional change (£* rne - ^observed/€true)- For 
objects with S/N = 10 and high intrinsic ICD (£true > 
6 %), the bias is small, with relative changes of < 34%, 
and the bias decreases for objects with increasing ICD. 
The bias increases for objects with lower intrinsic ICD, 
such that the relative change is 100 % for objects with 
S/N = 10 and £true = 2%. This is acceptable because 
the ICD increases by only a factor of 2 . Therefore, to our 
S/N = 10 limit, we can be confident that objects with 
observed high ICD have high intrinsic ICD, and objects 


with low observed ICD have low intrinsic ICD. 

In principle this S/N requirement applies to both bands 
used in Equation 2. However, for this particular study, 
the redder band (namely Hi 6 o) S/N is much higher in the 
majority of galaxies analyzed, such that this requirement, 
in practice, only applies to the bluer imaging band. It is 
also important to note that this S/N does not correspond 
to a single bluer band magnitude limit, because at fixed 
magnitude a larger object (with lower surface brightness) 
will have a lower S/N. For galaxies with an intrinsic ICD 
very near zero, and low surface brightness in the bluer 
band, it is possible for the ICD calculation to produce 
negative results when the background contributes signif¬ 
icantly to the fluctuations. Galaxies with an intrinsically 
larger ICD are more resistant to surface brightness dim¬ 
ming as the fluctuations of the galaxy dominate those of 
the background. 

After applying the S/N = 10 restriction in U 75 , we have 
a final, primary sample of 2601 galaxies. Unfortunately, 
this S/N requirement leads to a selection effect where we 
prefer galaxies with some level of active star-formation, 
as passively evolving galaxies often do not have enough 
blue light to be included in this analysis. The effect of 
this selection bias is discussed further in Section 4.1. 

3.3. Error Associated with ICD 

Papovich et al. (2003) defines the statistical error asso¬ 
ciated with the ICD, assuming that the background pix¬ 
els (those not associated with any galaxy) are normally 
distributed, as 


1ST • 

1 y pxx 

m = jr. --• ( 3 ) 

iV pxx pxx 

E (hi -P) 2 - E (B2,i - a£?i,i) 2 

i=1 i= 1 

Papovich et al. (2005) note that Equation 3 underesti¬ 
mates the uncertainties associated with the ICD. They 
estimated that the uncertainties for Equation 2 in their 
data required corrections of factors of 6-7 to account for 
correlated noise in the images. Bond et al. (2011) find 
that even after the correction suggested by Papovich 
et al. (2005), the uncertainty on the ICD obtained for 
bright objects (those with low errors) are underestimated 
by as much as a factor of three. The dominant source 
of error being contributions from non-uniform sky back¬ 
grounds and PSF mismatch between the two passbands. 

For this work, we use the observed ICD values for dif¬ 
ferent levels of intrinsic ICD and galaxy S/N to estimate 
the error associated with each measurement. To estimate 
the systematic uncertainty we compare the median ICD 
fractional change (the dashed line in Figure 1) to the 
true value, and to estimate the statistical uncertainty we 
compute the median absolute deviation. 

In all cases, galaxies with very high (i 775 S/N > 100) 
S/N have low (< 1 %) systematic uncertainties and very 
low (< 0.2%) statistical uncertainties. The closer to 
our U 75 S/N cutoff (U 75 S/N = 10) the systematic error, 
described in the previous section, increasingly dominates. 
The statistical uncertainty remains very low ^ 0.2% and 
so we neglect any effect due to the statistical variation 
in the ICD. 
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4. THE RELATION BETWEEN GALAXY 
PROPERTIES AND INTERNAL COLOR 
DISPERSION 

For the majority of our results we choose to use 
£( 2775 , H 160 ) as our primary diagnostic. For the galaxies 
in our sample, these two filters span the Balmer/4000 
A break over z ~ 1.5 — 3, and allow us to compare 
the rest-frame UV (traced by the 2775 band) to the rest- 
frame optical (traced by the H 160 band). This provides 
us the ability to investigate recent star formation (with 
2775 ) and past-averaged star formation (with Hi 60 ). In 
several sections (Section 4.1 and Section 5.4) we make 
use of different versions of the ICD, namely £(V 606 , J125) 
which allows us to compare slightly different parts of 
each galaxy’s spectral energy distribution to better un¬ 
derstand the results given by £( 2775 , Hi 6 o)- 

The median ICD for the full sample is £( 2775 , Hi 60 ) ~ 
5%. In the sections that follow we refer to galaxies with 
ICD value above and below this median as high and low 
ICD, respectively 

4.1. The Internal Color Dispersion vs. Galaxy Mass 

Stellar mass is one of the most fundamental properties 
of galaxies. The stellar mass of a galaxy correlates with 
many of its physical properties including overall color 
(Griitzbauch et al. 2011) and SFR (e.g., Daddi et al. 
2007; Magdis et al. 2010) to 2 ~ 3 (and perhaps beyond 
e.g., Salmon et al. 2014). As such, it is possible to gain 
much insight into the formation and evolution of galaxies 
through the rates of change of the ICD as functions of 
stellar mass. 

Figure 2 shows the £(z775,Hi 60 ) as a function of stel¬ 
lar mass for the galaxies in our sample. We find the 
median £( 2775 , H 160 ) values in the three low stellar mass 
bins are between 3.5% — 4.5%, the median ICD value in 
the moderate stellar mass bin is above 8%, and at stel¬ 
lar masses Log M/M 0 > 11, we find that galaxies have 
a median ICD value about 9%. If we take into account 
the systematic underestimations of low ICD values, dis¬ 
cussed in Section 3.3, we find little relative change in the 
median £(2775, Hieo) between bins, in that the moderate 
stellar mass bins still show higher £( 2775 , H 160 ) medians 
than the bins of lower stellar masses. 

To illustrate the ICD, Figure 3 shows galaxy images in 
the £( 2775 , Hi 6 o)-mass plane. The galaxies are selected 
at random from all of the galaxies contained in each bin 
of mass and £( 2775 , Hieo)- The RGB images are combi¬ 
nations of H 160 Ji 25 , and Z 775 . Galaxies with elevated 
ICD values (top rows) and in the moderate stellar mass 
range appear significantly different than the galaxies in 
other stellar mass bin or galaxies with lower ICD values. 
Visually, galaxies with high ICD have morphologies with 
well-separated regions of rest-frame UV (blue) and rest- 
frame optical (red) light. Galaxies with low ICD values 
tend to have smooth morphologies and uniform colors, 
suggesting fewer regions of separated red and blue emis¬ 
sion. 

4.1.1. High Mass Galaxies 

The highest stellar mass galaxies (Log M/M 0 >11) 
deserve special attention. Because the S/N threshold ap¬ 
plies only to the bluer band used in the ICD calculation, 
galaxies which are well detected in H 160 may not be well 


detected in 2775 . Specifically, massive galaxies often have 
the reddest 2775 — H 160 colors, so we expect that many of 
them have high Hi 60 S/N and low 2775 S/N. The concern 
is that these objects could have higher intrinsic ICD, but 
are excluded because of our signal to noise requirement 
(see Section 3.2). 

To investigate this possibility, we compare the mean 
2775 — Hi 60 color of the high stellar mass galaxies in¬ 
cluded and excluded from our sample. Galaxies which 
have S/N ^ 775 > 10 have (2775 — Hieo) = 2.6 mag while 
galaxies with S/N ^ 775 < 10 have (2775 — H 160 ) = 3.4 mag. 
This suggests that the excluded high stellar mass galax¬ 
ies simply do not have enough blue light to produce high 
ICD values. Papovich et al. (2005) simulate just such 
a situation where the authors find that extremely red 
galaxies, even with the addition of highly star-forming 
regions produce nearly no ICD as the older stellar pop¬ 
ulations simply overwhelm any diversity. For the seven 
high stellar mass galaxies with ^775 — Hieo > 3 mag in our 
sample, four have high £( 2775 , H 160 ); however, a visual 
inspection provides additional insight into these objects 
(discussed below). 

A visual inspection of the 26 most massive galaxies, 
shown in Figure 4 shows that the majority of these ob¬ 
jects have smooth light profiles in both U 75 , and Hieo- 
There are several exceptions which fall into two broad 
categories: galaxies which appear to have faint, nearby 
objects detected in 2775 but not Hi 60 (e.g., 16375, 17205 
and 16934). Because H 160 serves as the object detec¬ 
tion image, it is possible that these faint objects lie 
within the Hieo pertrosian radius used for the ICD cal¬ 
culation, contributing significantly to the measured ICD. 
While the origin of these nearby objects is unclear, pos¬ 
sible scenarios include a chance alignment or a disrupted 
satellite. And, galaxies in the highest-mass subsample 
with extended structures, showing distinct red/blue re¬ 
gions and likely individual clumps of star-formation em¬ 
bedded within the galactic structures (e.g., 282, 15639 
and 15261). In both cases, these features contrast with 
the Hi 60 dominated central region producing the large 
£(i 775 ,Hi 6 o) value. The blue central regions in the 
2775 — H 160 color maps reflects small mismatches in the 
PSF matching (see van der Wei et al. 2012 for a descrip¬ 
tion) and do not contribute significantly to the ICD as 
they are small compared to the total flux of the galaxy. 

After a visual inspection, if we remove the galaxies 
with neighboring objects from the high stellar mass bin 
the median £( 2775 , H 160 ) becomes ~ 7%. However, if we 
instead place a more stringent i*? 75 S/N requirement (U 75 
S/N > 30) on the high stellar mass galaxies, we find the 
median £( 2775 , H 160 ) further drops to 3.3%. 

Because it is unclear whether or not the faint objects 
neighboring the primary galaxy are physically associated 
with the primary, or are chance alignments along the 
line of sight, we do not exclude these objects and do 
not adjust the 2775 S/N requirement. Instead, we note 
that these are objects are infrequent in the sample, and 
so do not affect our conclusions. So while the median 
£( 2775 ,Hieo) in the high stellar mass bin may appear to 
be very similar to the moderate stellar mass bins, it is 
likely that the median is being inflated due to inclusion 
of galaxies with neighboring objects or galaxies with low 
i 77 5 S/N and high Hi 60 S/N. 
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Log Mass (M©) Log Mass (M©) 

Data Q Quartiles — Medians 


Fig. 2.— The ICD as a function of galaxy stellar mass. Individual galaxies are shown as gray points. The galaxies are then binned 
together by stellar mass in bins of 0.5 dex (except for the highest mass bin which is 1.0 dex) beginning with Log M/M 0 = 8.5. The upper 
and lower quartiles of ICD are indicated by the blue boxes. The median of each bin is indicated by the red line. The upper and lower black 
bars show the 95th and 5th percentiles. The ICD for galaxies in our sample is low for low stellar masses, peaks in the intermediate stellar 
masses, and returns to low values at high stellar masses. 


4.2. The Internal Color Dispersion vs. Sersic Profile 

Many previous works (e.g., Blanton et al. 2003; Bell 
et al. 2012; van der Wei et al. 2012) use the Sersic index 
to classify a galaxy’s morphology. The Sersic index pa¬ 
rameterizes a galaxy’s surface brightness profile, where 
n = 4 corresponds to a classical bulge, r 1 / 4 -law profile 
(de Vaucouleurs 1948) and n = 1 corresponds to the ex¬ 
ponential profile of a galactic disk. We refer to galaxies 
with n < 2 as disk-dominated, and those with n > 2 as 
bulge-dominated. 

Figure 5 shows the f ( 2775 , H 160 ) as a function of stel¬ 
lar mass, separated based on their Sersic indices. Again, 
we divide the sample into mass bins, where each bin is 
0.5 dex wide or larger to ensure there are at least 10 
objects in each bin. Galaxies with high ICD prefer a 
lower Sersic index (disk dominated). Inside of the in¬ 
dividual Sersic bins, galaxies with intermediate stellar 
mass (10 < Log M/M 0 < 11) display the largest ICD, 
but as the Sersic index increases, the upper envelope of 
ICD values decreases: bulge-dominated (n > 2) galaxies 
have low ICD. For disk-dominated galaxies (n < 2) the 
median galaxy across all masses has a high ICD value. 

The systematic underestimation described in Sec¬ 
tion 3.3, has very little effect on the preference of galax¬ 
ies with low Sersic index and intermediate stellar mass 
to have high £( 2775 , Hi6o)- This is because there is very 
little relative change in the medians of the bins due to 
the underestimation of the ICD. 


There is a well known trend between the Sersic index 
and stellar mass, with the Sersic index increasing with in¬ 
creasing stellar mass, This relationship is apparent both 
at low (e.g., Dutton 2009) and high (e.g., Franx et al. 
2008; Bell et al. 2012) redshift. We find a similar trend. 
The three low Sersic bins lack (or have very few) galaxies 
with masses greater than Log M/M 0 ~ 11. In fact, it is 
only the highest Sersic bin which shows any significant 
number high stellar mass galaxies. 

4.3. The Internal Color Dispersion vs. Star Formation 

Rate 

The ICD is sensitive to heterogeneities in the stellar 
populations of galaxies. One possible cause of these het¬ 
erogeneities is unobscured, young stellar populations sep¬ 
arated from older, passively evolving populations. Be¬ 
cause young populations are needed to produce measur¬ 
able ICD, the star formation rate (SFR) of a galaxy could 
be closely related to the ICD. Figure 6 shows the specific 
star formation rate (sSFR), the star formation rate nor¬ 
malized by the galaxy’s stellar mass, as a function of 
£(^775, H160). 

We find that the £( 2775 , H 160 ) of a galaxy in not di¬ 
rectly tied to its sSFR. Galaxies with low f ( 2775 , Hieo) 
values show no correlation with sSFR, however, galaxies 
with high £(i 775 ,Hi 6 o) values show predominately “nor¬ 
mal” sSFRs. For galaxies with low sSFRs we expect 
low £(i 775 ,Hi 6 o) values as such galaxies have few young 
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Fig. 3.— Galaxy montage showing examples of galaxies as a function of their location on the ICD - stellar-mass plane. Columns are 
bins of stellar mass while rows are bins of ICD. The RGB images are combinations of Hi6o, J 125 , and 2775 . Galaxies with elevated ICD 
values (£( 1775 , H 160 ) > 15%) and stellar masses Log M/Mq 10 — 11 appear significantly different than the galaxies in either other stellar 
mass bin or galaxies with depressed ICD values. Each cutout is 2.5" on the sky regardless of the redshift of the object, which corresponds 
to ~ 20 kpc over the redshift range. Columns are bins of stellar mass while rows are bins of ICD. Galaxies filling the grid are selected at 
random from all galaxies falling into the bin. Empty squares indicate regions in the £( 2775 , Hi 6 o)-mass space where there are no galaxies 
in our sample. The labels represent the centers of either the ICD or stellar mass bins. So a galaxy in the bottom left corner (e.g. 23594) 
could have a stellar mass 8.5 < Log M/Mq < 9 and and ICD of 0% < £( 2775 , H 160 ) < 10%. 


stellar populations producing the blue light needed to 
contrast the redder light from more evolved stellar popu¬ 
lations. Galaxies showing very high sSFRs, often due to 
re-radiation by dust in the IR, also have low £( 2775 , Hi 60 ) 
suggesting that dust variations play a small role in the 
production of ICD. 

Again, this result is robust even after accounting for 
the systematic underestimate of the low £( 2775 , H 160 ) 
galaxies. Shifting the lowest £( 2775 ,^ 60 ) galaxies up¬ 
ward to higher £( 2775 , Hi 60 ) values does not produce a 
trend with sSFR. 

5. POTENTIAL DRIVERS OF THE ICD 

As discussed previously, the ICD is sensitive to het¬ 
erogeneous stellar population variations that give rise to 
differences in the UV-optical flux ratio within a galaxy. 
In this section we discuss possible scenarios which could 
lead to high ICD. Specifically, we investigate whether 
merger activity, the galaxy’s physical size, intrinsic color 


gradients, effects due to redshift, or heterogeneous ex¬ 
tinction drive increased ICD. 


5.1. Mergers 

Based on morphological asymmetries, Papovich et al. 
(2005) suggested that some of the galaxies with high 
ICD are likely the result of merger activity. To inves¬ 
tigate this, we use non-parametric morphological indica¬ 
tors, specifically the Gini coefficient, G, and the second 
order moment of the brightest 20 % of a galaxy’s flux, 
M 20 , to better understand the role of galaxy mergers in 
producing high ICD. Many studies demonstrate these pa¬ 
rameter’s usefulness in quantifying galaxy morphology in 
large samples at both low and high redshifts (e.g., Lotz 
et al. 2004; Abraham et al. 2007; Law et al. 2007, 2012 ; 
Wang et al. 2012 ; Lee et al. 2013). 
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Fig. 4 .— Galaxy montage of the high mass (11 > Log M/M©) galaxies. Each triplet of panels shows either the 2775, H160 imaging or 
the Z775 — H160 color map, along with the ID number and associated £(1775, H160) value. Each cutout is 2 . 5 " on the sky regardless of the 
redshift of the object, which corresponds to ~ 20 kpc over the redshift range. The panels are ordered from lowest to highest £(2775, Hi6o)- 
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Fig. 5. — Galaxies binned into four bins of Sersic index. The galaxies are then binned by stellar mass in bins of 0.5 dex beginning with 
Log M/M© = 8 . 5 . The upper and lower quartiles are indicated by the blue boxes. The medians of each bin is indicated by the red line. 
The upper and lower bars show the 95 th and 5 th percentiles. We find that galaxies with high ICD (greater than 5 %) are dominated by 
low Sersic indices in general suggesting they are either disk galaxies or have high extinction. 


Lotz et al. (2008) adopt the following classification 
scheme based on G-M 20 : 

Mergers : G > —0.14M 2 o + 0.33 (4) 

E/SO/Sa : G < -0.14M 20 + 0.33 and G > 0.14M 20 + 0.80 

Sb - Ir : G < -0.14M 20 + 0.33 and G < 0.14M 20 + 0.80 

Figure 7 shows the G and M 2 q values for the galax¬ 
ies in our sample broken into large bins of ICD. As 
can be seen, mergers are not a dominant driver of the 
ICD. The majority of galaxies with high ICD values 
(5% — 20 %) prefer disk-dominated morphologies, consis¬ 


tent with our findings between the ICD and Seric index 
above. This suggests that mergers are not a significant 
cause of elevated ICD, except, possibly, at the highest 
levels (£(^ 775 ?Hieo) > 20 %) where 21 % of galaxies are 
classified as mergers. 

Given the small number of galaxies in the highest 
£( 2775 ,Hieo) bin, we use a series of Monte Carlo simu¬ 
lations to estimate how often we would find a similar 
distribution simply by chance. We assign each galaxy a 
morphology based on their location in the G — M 2 q space 
while allowing their £(^ 775 , H 160 ) to be changed. We then 
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Fig. 6. — Galaxies’ specific star formation rate (sSFR) as a func¬ 
tion of £(^ 775 , Hi 6 o)- The gray points are the actual data values 
and the red line is the running median value. Galaxies which have 
been clipped for clarity are shown as black arrows. The shaded 
region shows the width of the star-forming “main sequence” for a 
1 O 1 O M 0 stellar mass galaxy at z — 2.25 taken from Speagle et al. 
(2014). We find that galaxies with the highest £( 2775 , Hi6o) have 
very close to the median sSFR while galaxies with the highest and 
lowest sSFRs lack many galaxies with very high £( 1775 , H 160 ) val¬ 
ues. 
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Fig. 7.— Morphological classification {G — M 20 ) of objects in 
our sample. Based on this classification scheme, objects above 
the dashed, green line are mergers, below the dashed, green line 
but above the blue, dotted line are elliptical galaxies and be¬ 
low both lines are spiral galaxies. The majority of galaxies with 
£(z 775 ,Hi 6 o) > 5% are not classified as mergers. 


select galaxies with £( 2775 , H 160 ) > 20% and record the 
fraction classified as mergers. We find that galaxies in 
the highest ICD bin are at least 21% mergers only 1.32% 
of the time. Or, given just random chance, 98.7% of the 
time we find fewer mergers in the highest ICD bin than 
what is actually measured. 

5.2. Physical Size 


Fig. 8. — £(2775 , H 160 ) as a function of the galaxy’s effective radii 
(semi-major axis). Gray points are the individual galaxies. The red 
line is the median £( 2775 , H 160 ) value and the blue dashed represent 
the interquartile range, 75% and 25% respectively. The dashed 
vertical line shows the effective radius corresponding to the FWHM 
of the WFC3 PSF (FWHM = 0.2") at z = 3.5. We find no trend 
suggesting a preference of physical size for a galaxy with elevated 
7 %) £(^775, Hi 6 o). 


The ICD requires that we are able to resolve scales to 
a minimum of 0.5 kpc (Papovich et al. 2003). Below this, 
the resolution of the image suppresses any ICD that may 
be present. Here we also considered how the ICD values 
change with increasing physical size. 

Figure 8 shows £(^ 775 , Hieo) as a function of semi-major 
axis, R e . For galaxies with R e below 3 kpc the median 
ICD remains low 5%). Considering only the inter¬ 
mediate mass (10 < Log M/M 0 < 11 ) we find ICD val¬ 
ues remain low until ^2.5 kpc after which it increases. 
In either case, we find more galaxies with large ICD at 
large R e . However, this is most likely a selection effect 
where we are more able to detect ICD in extended ob¬ 
jects. Furthermore, placing the galaxies on a size-mass 
relation provides a similar result. At a fixed stellar mass, 
we find more elevated £( 2775 , H 160 ) galaxies with larger 
radii than at small radii. 

One of our findings is that the ICD values are low for 
lower stellar mass galaxies, Log M/M 0 < 10. Because 
these galaxies have lower R e compared to higher mass 
galaxies (e.g., Franx et al. 2008; Williams et al. 2010) we 
test if the trend between ICD and R e affects this finding. 
We remove all objects with R e < 2 kpc (at all stellar 
masses) and recompute how the median ICD changes 
as a function of mass. The median £( 2775 , H 160 ) for 
Log M/M 0 < 10 remains low, with £(i 7 75, Hi 60 )median — 
5%. The median £( 2775 , Hi 60 ) for moderate mass galax¬ 
ies remains elevated, £(i 7 75, H 160 ) me ^ an > 10%, for 
10 <C Log M/M 0 <C 11, and the £( 2775 , H 160 ) declines 
at the highest masses. Therefore, the relation between 
ICD and R e does not dominate our conclusions between 
ICD and stellar mass. 

5.3. Color Gradients 

Galaxies show a radial color gradient when the color of 
the light in an annulus on the outskirts differs from the 
color inside a central region (see Section 2.6). For galax¬ 
ies with large color gradients, this difference could po¬ 
tentially lead to large ICD values. Szomoru et al. (2011) 
found that the majority of a mass-selected sample from 
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Fig. 9.— The galaxy’s color gradient, A(z 775 — Hi 6 o), as a func¬ 
tion of £( 1775 , Hi 6 o)- Gray points are the data and the red line 
is the median color gradient as a function of ICD. We find that 
galaxies with high ICD have systematically lower color gradients, 
suggesting they have redder cores than outskirts. 



Fig. 10.— The fractional change between the £( 2775 , H 160 ) cal¬ 
culated with the center 0.18 // removed and the £( 2775 , H 160 ) cal¬ 
culated for the entire galaxy as a function of £( 1775 , H 160 ) for the 
entire galaxy. We find that as galaxies with their centers removed 
show an elevated amount of £( 2775 , Hi 6 o)- This is expected as the 
variations in the outskirts play a more significant role with the 
bright central region removed, and indicates the observed ICD is 
not due to an underlying color gradient. 


the HUDF show negative color gradients (redder cores 
and bluer outer regions). We find a similar result for the 
galaxies in our sample, with an average color gradient 
^(^775 — H 160 ) = — 0.12 and 66 % of the sample showing 
a negative color gradient. 

The relationship between color gradient and the ICD is 
complex. Figure 9 shows A (2775 — H 160 ) as a function of 
£( 2775 , Hieo)- We find that while galaxies with high ICD 
typically have negative color gradients, there is no corre¬ 
lation between the strength of the color gradient and the 
ICD. Therefore, while strong color gradients do not nec¬ 
essarily imply high ICD values, galaxies with high ICD 
appear to require larger color gradients, either positive 
or negative. 

To understand how much a color gradient contributes 
to the measured ICD, we perform the following test. We 
mask the central 0.18" 1.5 kpc) of all galaxies and 

recompute the ICD value. Figure 10 shows the relative 
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Fig. 11.— ICD of galaxies in our sample as a function of redshift. 
The galaxies are then binned together by redshift in bins of 0.5 
dex beginning with z — 1.5. The upper and lower quartiles are 
indicated by the blue boxes. The medians of each bin is indicated 
by the red line. For clarity, only a portion of the possible ICD 
range is shown. Top and bottom panel show £( 2775 , H 160 ) and 
£(V606 5 J 125 ) as a function of redshift, respectively. In both panels 
we do not see a significant increase in the median values of ICD. 
The inclusion of the bluer bands (V606 and J 125 ) acts as a probe 
for bandpass shifting effects. Because neither shows a significant 
trend, this suggests that the lack of increased values is not simply 
due to bandpass shifting effects. 

change in the ICD due to this masking. In the majority 
of galaxies with small £( 2775 , H 160 ) values, removing the 
nuclear region causes the ICD to decrease, while in galax¬ 
ies with larger ICD values there is little change. This is 
as expected as galaxies with small ICD values, masking 
out the central regions also removes some of the varia¬ 
tion. Therefore, while galaxies with strong color gradi¬ 
ents may have high ICD, this is not a requirement. The 
observed color gradients alone are insufficient to produce 
the ICD in the galaxies in our sample. 

5.4. Redshift Effects from Bandpass Shifting on the 

ICD 

The galaxies in our sample span a range of redshift 
from 1.5 < z < 3.5. As a result, the rest-frame light 
probed by both 2775 and H 160 varies between galaxies. 
For example, in galaxies at z = 1.5 — 2.5, Hieo probes 
rest-frame 4600 — 6400 A , while at higher redshifts, z > 
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Fig. 12. — The Internal Color Dispersion as a function of 
Z 775 — Hi 60 color. The colored boxes show the average attenua¬ 
tion value (Ay) of the points within the region. For the galaxies in 
our sample, we find no evidence of a correlation between the mea¬ 
sured stellar population diversity and the rest frame UV-optical 
color. We do find that galaxies with large ICD values also show 
large attenuation values. 

3, H 160 probes light shortward of 4000 A . Therefore, 
we have considered how the variation in rest-frame light 
could affect our ICD measurement. 

For galaxies at 2: < 2, we measure the ICD be¬ 
tween V 60 6 and J 125 , as these bandpasses observe sim¬ 
ilar rest-frame wavelengths as 2775 and H 160 for galax¬ 
ies at z > 2.2. Figure 11 shows ICD as a function of 
redshift where we have binned the galaxies by redshift. 
The top row of Figure 11 shows £( 2775 , H 160 ) and the 
bottom shows £(V 606 5 Ji2s)- Dor the comparable bins 
(z = 1.5 — 2 for £(V606? J 125 ) and 2 ~ 2.2 for £( 2775 , Hi 60 )) 
we find very little (< 0.1%) difference between the me¬ 
dian ICD values. This shows that as long as the imaging 
bands chosen to calculate the ICD span the rest frame 
Balmer/4000A break we calculate a similar ICD. 

For galaxies above 2 = 3 both 2775 and H 160 he blue- 
ward of the rest frame Balmer/4000 A break, and Bond 
et al. ( 2011 ) suggests that a sharp downturn would be ob¬ 
served as the ICD probes similar populations. We find no 
such downturn, suggesting that while the Balmer/4000 
A break is moving through H 160 , the observed trends in 
£( 2775 , H 160 ) as a function of stellar mass (Figure 2 ) is 
not simply due to the effect of redshift, but due to the 
individual characteristics of the galaxies. 

5.5. Effects of Dust Attenuation on the ICD 

As mentioned earlier, high ICD could result from spa¬ 
tial variations of stellar ages or dust attenuation. In this 
section, we investigate whether the ICD is related to the 
total dust content of our galaxies. 

If inhomogeneous attenuation (the reddening of selec¬ 
tive parts of a uniformly blue population) is the primary 
driver of the ICD then we would expect to find galaxies 
with both very red colors and high ICD values. A uni¬ 
formly blue galaxy would begin with a blue 2775 — H 160 
color and if we were to add a level of inhomogeneous dust 
then we would expect the ICD to rise, as there is now 
red light with which to contrast the blue and the overall 
color to shift toward the red. Similarly if a galaxy with 
a heterogeneous stellar populations has the young star 


Fig. 13. — The ratio of the total SFR to the SFR derived from 
observed UV luminosity, SFR-potal/SFRuv, as a function of ICD. 
The colored boxes show the average attenuation value (Ay) of the 
points within the region. Galaxies with IR dominated SFRs have 
smaller ICDs than those with a lower SFRxotai/SFRuv- At fixed 
SFR T otal/SFR uv , there is a large spread in derived attenuation. 
This suggests that a mix of stellar population ages drives the ICD 
and attenuation plays a minor role. 

forming regions attenuated then the ICD would go down 
as there is less blue light contrasting the red and the 
2775 — Hi 60 color again shifts toward the red. Finally, if a 
galaxy dominated by older stellar populations is dust at¬ 
tenuated then the red stellar populations become more 
red, shifting the 2775 — Hieo color, but the ICD would 
increase as there is now some contrast between stellar 
populations. 

In Figure 12 we show the ICD as a function of 2775 — 
H 160 color and attenuation value (Ay) as reported by 
FAST (see Section 2.4). There is little evidence for a 
correlations between either the ICD or Ay as we find 
high £( 2775 , H 160 ) galaxies for a wide range of colors and 
Ay. The ICD is low for galaxies with both high and low 
2775 — Hi 60 colors, but galaxies with moderate 2775 — Hi 60 
color (1 mag < 2775 — H 160 < 2 mag) show the largest 
range in the ICD. Therefore, it seems that galaxies with 
moderate ^ 775 —Hieo colors are required for high ICD, but 
color alone is not enough to predict high ICD. Papovich 
et al. (2005); Bond et al. (2011) argue this lack of direct 
correlation supports the idea that the ICD depends on 
spatial variations in the ages of the stellar populations 
and not dust attenuation. 

As a second investigation, we compare the ratio of the 
total galaxy SFR (SFR to tai ) 5 derived from the IR and 
UV emission (see Section 2.7), against the SFR inferred 
from the UV measurement (uncorrected for dust atten¬ 
uation). In this way, the ratio of SFR to tai/SFRuv is a 
measure of the amount of light from star formation that 
is attenuated by dust and reradiated in the IR. Figure 13 
shows there is no trend in SFR to tai/SFRuv as a function 
of £(^ 775 , Hi 6 o)- Galaxies with high £(^ 775 , H 160 ) span the 
entire range of SFR tota i/SFRuy. 

Taken together, this is evidence that spatial variations 
of dust is not the main driver of the ICD, and as such, 
the ICD requires spatial variations in stellar population 
ages in galaxies (perhaps with contributions from spatial 
variations in dust attenuation). This agrees with the 
conclusions from Papovich et al. (2003, 2005) and Bond 
et al. (2011). Both authors acknowledge that extinction 
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Fig. 14.— ICD as a function of number of clumps. The gray 
points are the individual galaxies. The red solid and blue dashed 
lines represent the average clump number in ICD bins and the 
average ICD value for different numbers of clumps. 


variations could play a role in providing variations in 
the blue starlight, but conclude that the effect is smaller 
than those produced by variations in segregated stellar 
populations of different mass-weighted ages. The ICD 
appears to driven mostly by variations in the ages of 
stellar populations in galaxies. 

5.6. Star Forming Clumps 

Papovich et al. (2005) suggest that in order to have 
ICD > 10%, young star-forming disk in galaxies are re¬ 
quired to have a surrounding older spheroid; for higher 
ICDs (> 20 — 25%), they suggest compact star-forming 
HII regions (discrete star-forming regions or clumps) are 
also required. This picture is corroborated in a recent 
study by Wuyts et al. ( 2012 ) who found that star-forming 
galaxies at z ~ 2, with stellar masses Log M/M© > 10 
have redder centers and bluer outer regions. 

Figure 14 shows £(^ 775 , Hi 60 ) as a function of the num¬ 
ber of clumps, using the catalog of Guo et al. (2014), 
see Section 2.8. The average £( 2775 , H 160 ) value steadily 
increases for increasing number of clumps but turns over 
after a galaxy contains at least four clumps. The average 
number of clumps for a given £( 2775 , H 160 ) value shows 
a similar increase with increasing £(i 775 ,Hi 60 ) but does 
not show a turnover. For galaxies with low (< 5%) ICD 
values, the mean number of clumps per galaxy is 1 . 2 ; 
for galaxies with higher ICDs (> 5%), the mean number 
of clumps is 2.1 Therefore, we find that if a galaxy has 
at least two clumps, on the average, it will have a high 
£( 2775 , Hi 60 ) value. This trend continues, where galaxies 
with four or more clumps show median ICD of greater 
than 10%. Therefore, the ICD is sensitive to the number 
of star-forming clumps. 

The clump catalog provides several clump detection 
levels (see Section 2.8). When we increase in the detec¬ 
tion threshold (UV luminosity > 5% instead of 1%) we 
find the mean number of clumps remains roughly con¬ 
stant (1 — 2) regardless of the £( 2775 , H 160 ) value. Also, 
while the total number of clumps per galaxy is reduced, 
we find the £(^ 775 , H 160 ) of galaxies with fewer than two 
clumps remains generally unchanged between the two 
UV luminosity thresholds. 

Because the star-forming clumps are expected to form 


in gas-rich disks (e.g. Keres et al. 2005; Dekel & Birnboim 
2006; Dekel et al. 2009b; Dekel & Burkert 2014), the ICD 
must trace galaxies which have a redder but still gaseous 
stellar disk. These individual clumps are thought to be 
relatively short lived (only 100 — 200 Myr; Wuyts et al. 
2012 ), which would correspond well with the lifetimes of 
hot, blue stars contributing to the UV light associated 
with high ICD. 

6 . ON BULGES, CLUMPS, AND THE ICD 

For the galaxies in this study, we find that as both the 
stellar mass and Sersic index increase, the ICD decreases 
(Figure 5). This is unsurprising as previous studies sug¬ 
gest galaxies with high Sersic indices have smooth light 
profiles with little color variation across them (e.g. van 
Dokkum et al. 2004; Daddi et al. 2004; Szomoru et al. 
2010 ) and these stellar populations produce little color 
dispersion as any variation is small relative to the total 
galaxy flux (Papovich et al. 2005). 

However, the galaxies with high ICD values have in¬ 
termediate stellar mass and are disk-dominated (n < 2 ). 
As mentioned previously, Papovich et al. (2005) found 
that for a galaxy to have high ICD a red bulge sur¬ 
rounded by a blue disk with HII regions (i.e. clumps) 
is required. Because high ICD is most prominent in 
disk-dominated morphologies and associated with star¬ 
forming clumps, the ICD is less sensitive to the red 
bulge component and more driven by the color varia¬ 
tion between the underlying stellar disk and highly star¬ 
forming clump. We see this in the previous comparison 
of the £(^ 775 ?Hi 6 o) computed with and without the cen¬ 
tral region (see Section 5.3), where we found that galax¬ 
ies with their nuclear region removed show an increased 
C(* 775 ,Hi 60 ) (Figure 10). 

Interpreted another way, the decrease of © 775 , Hieo) 
with increasing Sersic index suggests that there is a con¬ 
tinual suppression of the mechanism which causes high 
ICD (perhaps clumps) with spheroid formation. This 
shows that if galaxies are to form quiescent bulges or 
spheroids along with a star-forming disk or clumps (as is 
typical of classical Hubble Sequence galaxies) it is only 
happening at this moderate (10 < Log M/M© < 11) 
stellar mass range and with relatively small bulge con¬ 
tributions. Once a bulge contributes significantly to the 
galaxy’s morphology, the amount of ICD decreases dra¬ 
matically, suggesting that the bulge plays a role in the 
homogenization of the galaxy’s stellar populations. 

Two different effects could be associated with the for¬ 
mation of a bulge and subsequent reduction in the ob¬ 
served ICD. One theory is the “morphological quench¬ 
ing” (Martig et al. 2009) concept, by which the forma¬ 
tion of the a bulge stabilizes the gaseous disk against fur¬ 
ther fracturing and star formation. The second could be 
some feedback process associated with the bulge which 
removes any remaining gas from the disk, preventing new 
clump formation. 

The first step in the quenching process is to form a 
strong bulge component. One scenario for this process is 
for clumps to migrate through dynamical friction to co¬ 
alesce in the nuclei, and merge to form a bulge (e.g., 
Elmegreen et al. 2008; Dekel et al. 2009b, a; Ceverino 
et al. 2010). The galaxies with high ICD values also 
show negative color gradients (Figure 9) suggesting the 
central stars are older (or perhaps more attenuated by 
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dust) than the outer star-forming regions. Other studies 
(e.g., Guo et al. 2012; Wuyts et al. 2012, 2013; Adamo 
et al. 2013) find that central clumps are both older (by 
500 Myrs) and more attenuated than the clumps in 
the outskirts suggesting that migration plays a significant 
role. 

Of course, not all galaxies in the moderate stellar mass 
range have high ICD values. Either only some galaxies 
have periods of high ICD or all galaxies go through such 
a phase, but the period of elevated ICD is very short 
lived. If all galaxies experience a period of elevated ICD 
at some point during their evolution, then the character¬ 
istic time scale of high ICD is the product of the time 
spanned by the look-back time across the redshift range 
(2.4 Grys) and the fraction of sources with high ICD. 
In this work, 21% of galaxies have £( 2775 , Hi6o) > 10%. 
Therefore, if all galaxies go through a high ICD phase, 
then the timescale of high ICD is approximately 375 
Myrs. This timescale is similar to the predicted clump 
survival timescale of 200 — 250 Myrs (Wuyts et al. 2012). 

Regardless of the mechanism at work, there is a dis¬ 
tinct suppression in the ICD of galaxies once a bulge 
is formed. Galaxies which contain highly star forming 
clumps match well the description of these high ICD ob¬ 
jects and clump migration could be a primary cause of 
ICD suppression. It will take further investigation in the 
future to better associate galaxies showing high ICD with 
star forming clumps. 

7. SUMMARY 

In this work, we investigate the dependence of the rest- 
frame UV and optical ICD on other galaxy properties. 
We define the ICD as the ratio of the square difference of 
the image flux - intensity values about the mean galaxy 
color to the square of the total image flux. It is a measure 
of the amount of variation in the stellar diversity of a 
galaxy; that is separated red and blue stellar populations. 
Using imaging from CANDELS of GOODS-South and 
the HUDF, we select galaxies between 1.5 < z < 3.5, a 
Hieo magnitude brighter than 25 mag and a requirement 
that the signal-to-noise be sufficient enough to provide 
an accurate measurement of the ICD (see Section 2.3). 
This selection yields 2601 galaxies. 

We present our results from comparing the ICD to the 
stellar mass of the galaxies (Section 4.1), the morphol¬ 
ogy based on Serisc index (Section 4.2) and star for¬ 
mation rates (Section 4.3). We then investigate possi¬ 
ble causes of the ICD including: mergers (Section 5.1), 
physical size (Section 5.2), color gradients (Section 5.3), 
band-pass shifting effects (Section 5.4), attenuation (Sec¬ 
tion 5.5) and clumpy star forming regions (Section 5.6). 
We attribute ICD to heterogeneous stellar populations 
(possibly clumps) with a possible contribution from dust 
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attenuation. 

The major results of this work, with accompanying 
interpretations, are as follows: 

• Galaxies have the largest amount of stellar popu¬ 
lation diversity, as measured by the ICD, at stellar 
masses 10 < Log M/M 0 < 11. This stellar mass 
range corresponds well to the range in which local 
galaxies experience a decrease in their star forma¬ 
tion efficiency, and galaxies at z ^ 2 which are 
found to show star forming clumps. Therefore, the 
stellar population diversity driving the ICD could 
be discrete, highly star-forming regions contrasted 
against an older, redder, underlying population. 

• High ICD is preferentially found in galaxies with 
extended light profiles (Serisc index < 2). This 
suggests that galaxies only have diverse stellar pop¬ 
ulations when there is little or no bulge present. 

• There is a continual suppression in the mechanism 
which causes high ICD with increasing Serisc in¬ 
dex. This could be due formation of clumps, which 
initially increase the ICD, and then migrate to¬ 
ward the center and coalescence into a quiescent 
bulge. The presence of this bulge then stabilizes 
the disk again fragmentation and new clump for¬ 
mation, which suppresses the ICD. 

The results suggests that some fraction of galaxies 
with high ICD are a result of the combination of struc¬ 
tures such as quiescent bulges and star-forming (perhaps 
clump) disks in galaxies. If this corresponds to the for¬ 
mation of the spiral galaxies indicative of the Hubble 
Sequence, then this phase occurs only in galaxies with 
stellar masses 10 < Log M/M 0 <11 and when the bulge 
component remains relatively small (n < 2). 
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